Cysteine hyperoxidation rewires communication pathways in the nucleosome and destabilizes the dyad.

Gene activity is tightly controlled by reversible chemical modifications called epigenetic marks, which are of various types and modulate gene accessibility without affecting the DNA sequence. Despite an increasing body of evidence demonstrating the role of oxidative-type modifications of histones in gene expression regulation, there remains a complete absence of structural data at the atomistic level to understand the molecular mechanisms behind their regulatory action. Owing to µs time-scale MD simulations and protein communication networks analysis, we describe the impact of histone H3 hyperoxidation (i.e., S-sulfonylation) on the nucleosome core particle dynamics. Our results reveal the atomic-scale details of the intrinsic structural networks within the canonical histone core and their perturbation by hyperoxidation of the histone H3 C110. We show that this modification involves local rearrangements of the communication networks and destabilizes the dyad, and that one modification is enough to induce a maximal structural signature. Our results suggest that cysteine hyperoxidation in the nucleosome core particle might favor its disassembly.

Design and validation of algorithms to identify venous thromboembolism in the French National Healthcare Database.

PURPOSE: This paper aims to introduce an algorithm designed to identify Venous Thromboembolism (VTE) in the French National Healthcare Database (SNDS) and to estimate its positive predictive value. METHODS: A case-identifying algorithm was designed using SNDS inpatient and outpatient encounters, including hospital stays with discharge diagnoses, imaging procedures and drugs dispensed, of French patients aged at least 18 years old to whom baricitinib or Tumor Necrosis Factor Inhibitors (TNFi) were dispensed between September 1, 2017, and December 31, 2018. An intra-database validation study was then conducted, drawing 150 cases identified as VTE by the algorithm and requesting four vascular specialists to assess them. Patient profiles used to conduct the case adjudication were reconstituted from de-identified pooled and formatted SNDS data (i.e., reconstituted electronic health records-rEHR) with a 6-month look-back period prior to the supposed VTE onset and a 12-month follow-up period after. The positive predictive value (PPV) with its 95% confidence interval (95% CI) was calculated as the number of expert-confirmed VTE divided by the number of algorithm-identified VTE. The PPV and its 95% CI were then recomputed among the same patient set initially drawn, once the VTE-identifying algorithm was updated based on expert recommendation. RESULTS: For the 150 patients identified with the first VTE-identifying algorithm, the adjudication committee confirmed 92 cases, resulting in a PPV of 61% (95% CI = [54-69]). The final VTE-identifying algorithm including expert suggestions showed a PPV of 92% (95% CI = [86-98]) with a total of 87 algorithm-identified cases, including 80 retrieved from the 92 confirmed by experts. CONCLUSION: The identification of VTE in the SNDS is possible with a good PPV.

Binding to nucleosome poises human SIRT6 for histone H3 deacetylation.

Sirtuin 6 (SIRT6) is an NAD+-dependent histone H3 deacetylase that is prominently found associated with chromatin, attenuates transcriptionally active promoters and regulates DNA repair, metabolic homeostasis and lifespan. Unlike other sirtuins, it has low affinity to free histone tails but demonstrates strong binding to nucleosomes. It is poorly understood how SIRT6 docking on nucleosomes stimulates its histone deacetylation activity. Here, we present the structure of human SIRT6 bound to a nucleosome determined by cryogenic electron microscopy. The zinc finger domain of SIRT6 associates tightly with the acidic patch of the nucleosome through multiple arginine anchors. The Rossmann fold domain binds to the terminus of the looser DNA half of the nucleosome, detaching two turns of the DNA from the histone octamer and placing the NAD+ binding pocket close to the DNA exit site. This domain shows flexibility with respect to the fixed zinc finger and moves with, but also relative to, the unwrapped DNA terminus. We apply molecular dynamics simulations of the histone tails in the nucleosome to show that in this mode of interaction, the active site of SIRT6 is perfectly poised to catalyze deacetylation of the H3 histone tail and that the partial unwrapping of the DNA allows even lysines close to the H3 core to reach the enzyme.

Synthesis of Original Cyclic Dinucleotide Analogues Using the Sulfo-click Reaction.

The stimulator of interferon genes (STING) protein plays a crucial role in the activation of the innate immune response. Activation of STING is initiated by cyclic dinucleotides (CDNs) which prompted the community to synthesize structural analogues to enhance their biological properties. We present here the synthesis and biological evaluation of four novel CDN analogues composed of an N-acylsulfonamide linkage. These CDNs were obtained in high overall yields via the sulfo-click reaction as a key step.

Resolving a guanine-quadruplex structure in the SARS-CoV-2 genome through circular dichroism and multiscale molecular modeling.

The genome of SARS-CoV-2 coronavirus is made up of a single-stranded RNA fragment that can assume a specific secondary structure, whose stability can influence the virus's ability to reproduce. Recent studies have identified putative guanine quadruplex sequences in SARS-CoV-2 genome fragments that are involved in coding for both structural and non-structural proteins. In this contribution, we focus on a specific G-rich sequence referred to as RG-2, which codes for the non-structural protein 10 (Nsp10) and assumes a guanine-quadruplex (G4) arrangement. We provide the secondary structure of RG-2 G4 at atomistic resolution by molecular modeling and simulation, validated by the superposition of experimental and calculated electronic circular dichroism spectra. Through both experimental and simulation approaches, we have demonstrated that pyridostatin (PDS), a widely recognized G4 binder, can bind to and stabilize RG-2 G4 more strongly than RG-1, another G4 forming sequence that was previously proposed as a potential target for antiviral drug candidates. Overall, this study highlights RG-2 as a valuable target to inhibit the translation and replication of SARS-CoV-2, paving the way towards original therapeutic approaches against emerging RNA viruses.

Robust AMBER Force Field Parameters for Glutathionylated Cysteines.

S-glutathionylation is an oxidative post-translational modification, which is involved in the regulation of many cell signaling pathways. Increasing amounts of studies show that it is crucial in cell homeostasis and deregulated in several pathologies. However, the effect of S-glutathionylation on proteins' structure and activity is poorly understood, and a drastic lack of structural information at the atomic scale remains. Studies based on the use of molecular dynamics simulations, which can provide important information about modification-induced modulation of proteins' structure and function, are also sparse, and there is no benchmarked force field parameters for this modified cysteine. In this contribution, we provide robust AMBER parameters for S-glutathionylation, which we tested extensively against experimental data through a total of 33 µs molecular dynamics simulations. We show that our parameter set efficiently describes the global and local structural properties of S-glutathionylated proteins. These data provide the community with an important tool to foster new investigations into the effect of S-glutathionylation on protein dynamics and function, in a common effort to unravel the structural mechanisms underlying its critical role in cellular processes.

Impact of Cardiovascular Comorbidities on the Effectiveness and Safety of Bevacizumab in Older Patients with Metastatic Colorectal Cancer.

BACKGROUND: Cardiovascular comorbidities are not contraindications of bevacizumab for metastatic colorectal cancer. OBJECTIVE: We aimed to evaluate the impact of cardiovascular comorbidities before bevacizumab treatment on overall survival and cardiovascular safety in older patients with metastatic colorectal cancer. METHODS: A 2009-2015 cohort of patients with metastatic colorectal cancer aged ≥ 65 years administered first-line bevacizumab was extracted from the French healthcare reimbursement claims database. Baseline heart failure, hypertension, and venous/arterial thromboembolic events were identified. The 36-month overall survival rate was evaluated using the Kaplan-Meier method, and the impact of cardiovascular comorbidities on the 36-month overall survival using a time-dependent, multivariable, Cox proportional hazards model. The 36-month cumulative incidence of cardiovascular events, and the impact of cardiovascular comorbidities on the likelihood of cardiovascular events were evaluated using the Fine and Gray model, with death as a competing risk. RESULTS: We included 9222 patients (56.4% male; median age 73 years). Two-thirds (66.7%) had baseline cardiovascular comorbidities. The median 36-month overall survival was 20.4 [95% confidence interval (CI) 19.9-21.0] and 21.8 [95% CI 21.1-22.6] months in patients with and without cardiovascular comorbidities, respectively. Age ≥ 75 years, dependency in activities of daily living, radiotherapy, and another targeted therapy were identified as death risk factors, but not cardiovascular comorbidities. At 36 months, cardiovascular events had occurred in 60.2% [95% CI 58.9-61.4] and 44.1% [95% CI 42.3-45.9] of patients with and without cardiovascular comorbidities. Baseline venous thrombosis, female, three or more cardiovascular medications, another targeted therapy, and more than six bevacizumab injections were identified as risk factors for cardiovascular events. CONCLUSIONS: In clinical practice, cardiovascular comorbidities before administering bevacizumab to older patients with metastatic colorectal cancer impacted the cardiovascular safety, but not overall survival. Unless they limit functional independency, older patients with cardiovascular comorbidities should be treated with bevacizumab under close monitoring.

Mechanism of the Covalent Inhibition of Human Transmembrane Protease Serine 2 as an Original Antiviral Strategy.

The Transmembrane Protease Serine 2 (TMPRSS2) is a human enzyme which is involved in the maturation and post-translation of different proteins. In addition to being overexpressed in cancer cells, TMPRSS2 plays a further fundamental role in favoring viral infections by allowing the fusion of the virus envelope with the cellular membrane, notably in SARS-CoV-2. In this contribution, we resort to multiscale molecular modeling to unravel the structural and dynamical features of TMPRSS2 and its interaction with a model lipid bilayer. Furthermore, we shed light on the mechanism of action of a potential inhibitor (nafamostat), determining the free-energy profile associated with the inhibition reaction and showing the facile poisoning of the enzyme. Our study, while providing the first atomistically resolved mechanism of TMPRSS2 inhibition, is also fundamental in furnishing a solid framework for further rational design targeting transmembrane proteases in a host-directed antiviral strategy.

Predicting the Three-Dimensional Structure of the c-KIT Proto-Oncogene Promoter and the Dynamics of Its Strongly Coupled Guanine Quadruplexes.

Guanine quadruplexes (G4s) play essential protective and regulatory roles within cells, influencing gene expression. In several gene-promoter regions, multiple G4-forming sequences are in close proximity and may form three-dimensional arrangements. We analyze the interplay among the three neighboring G4s in the c-KIT proto-oncogene promoter (WK1, WSP, and WK2). We highlight that the three G4s are structurally linked and their cross-talk favors the formation of a parallel structure for WSP. Relying on all-atom molecular dynamic simulations exceeding the µs time scale and using enhanced sampling methods, we provide the first computationally resolved structure of a well-organized G4 cluster in the promoter of a crucial gene involved in cancer development. Our results indicate that neighboring G4s influence their mutual three-dimensional arrangement and provide a powerful tool to predict and interpret complex DNA structures that can ultimately be used as a starting point for drug discovery.

How SARS-CoV-2 Alters the Regulation of Gene Expression in Infected Cells.

Nonstructural accessory proteins in viruses play a key role in hijacking the basic cellular mechanisms, which is essential to promote the virus survival and evasion of the immune system. The immonuglobulin-like open reading frame 8 (ORF8) protein expressed by SARS-CoV-2 accumulates in the nucleus and may influence the regulation of the gene expression in infected cells. In this contribution, by using microsecond time-scale all-atom molecular dynamics simulations, we unravel the structural bases behind the epigenetic action of ORF8. In particular, we highlight how the protein is able to form stable aggregates with DNA through a histone tail-like motif, and how this interaction is influenced by post-translational modifications, such as acetylation and methylation, which are known epigenetic markers in histones. Our work not only clarifies the molecular mechanisms behind the perturbation of the epigenetic regulation caused by the viral infection but also offers an unusual perspective which may foster the development of original antivirals.

Molecular Basis of the pH-Controlled Maturation of the Tick-Borne Encephalitis Flavivirus.

Flaviviruses are enveloped viruses causing high public concerns. Their maturation spans several cellular compartments having different pH. Thus, complex control mechanisms are in place to avoid premature maturation. Here we report the dynamical behavior at neutral and acidic pH of the precursor of the membrane fusion protein E of tick-borne encephalitis, showing the different stabilizations of the E dimer and the role played by the small fusion-assisting protomer (pr). The comprehension, at atomic resolution, of the fine regulation of viral maturation will be fundamental to the development of efficient strategies against emerging viral threats.

Revealing the Molecular Interactions between Human ACE2 and the Receptor Binding Domain of the SARS-CoV-2 Wild-Type, Alpha and Delta Variants.

After a sudden and first spread of the pandemic caused by the novel SARS-CoV-2 (Severe Acute Respiratory Syndrome-Coronavirus 2) wild-type strain, mutants have emerged which have been associated with increased infectivity, inducing surges in the contagions. The first of the so-called variants of concerns, was firstly isolated in the United Kingdom and later renamed Alpha variant. Afterwards, in the middle of 2021, a new variant appeared called Delta. The latter is characterized by the presence of point mutations in the Spike protein of SARS-CoV-2, especially in the Receptor Binding Domain (RBD). When in its active conformation, the RBD can interact with the human receptor Angiotensin-Converting Enzyme 2 (ACE2) to allow the entry of the virions into cells. In this contribution, by using extended all-atom molecular dynamic simulations, complemented with machine learning post-processing, we analyze the changes in the molecular interaction network induced by these different strains in comparison with the wild-type. On one hand, although relevant variations are evidenced, only limited changes in the global stability indicators and in the flexibility profiles have been observed. On the other hand, key differences were obtained by tracking hydrophilic and hydrophobic molecular interactions, concerning both positioning at the ACE2/RBD interface and formation/disruption dynamic behavior.

Understanding the Interactions of Guanine Quadruplexes with Peptides as Novel Strategies for Diagnosis or Tuning Biological Functions.

Guanine quadruplexes (G4s) are nucleic acid structures exhibiting a complex structural behavior and exerting crucial biological functions in both cells and viruses. The specific interactions of peptides with G4s, as well as an understanding of the factors driving the specific recognition are important for the rational design of both therapeutic and diagnostic agents. In this review, we examine the most important studies dealing with the interactions between G4s and peptides, highlighting the strengths and limitations of current analytic approaches. We also show how the combined use of high-level molecular simulation techniques and experimental spectroscopy is the best avenue to design specifically tuned and selective peptides, thus leading to the control of important biological functions.

Modeling the Enzymatic Mechanism of the SARS-CoV-2 RNA-Dependent RNA Polymerase by DFT/MM-MD: An Unusual Active Site Leading to High Replication Rates.

Viral infection relies on the hijacking of cellular machineries to enforce the reproduction of the infecting virus and its subsequent diffusion. In this context, the replication of the viral genome is a key step performed by specific enzymes, i.e., polymerases. The replication of SARS-CoV-2, the causative agent of the COVID-19 pandemics, is based on the duplication of its RNA genome, an action performed by the viral RNA-dependent RNA polymerase. In this contribution, by using highly demanding DFT/MM-MD computations coupled to 2D-umbrella sampling techniques, we have determined the chemical mechanisms leading to the inclusion of a nucleotide in the nascent viral RNA strand. These results highlight the high efficiency of the polymerase, which lowers the activation free energy to less than 10 kcal/mol. Furthermore, the SARS-CoV-2 polymerase active site is slightly different from those usually found in other similar enzymes, and in particular, it lacks the possibility to enforce a proton shuttle via a nearby histidine. Our simulations show that this absence is partially compensated by lysine whose proton assists the reaction, opening up an alternative, but highly efficient, reactive channel. Our results present the first mechanistic resolution of SARS-CoV-2 genome replication at the DFT/MM-MD level and shed light on its unusual enzymatic reactivity paving the way for the future rational design of antivirals targeting emerging RNA viruses.

Abiraterone acetate versus docetaxel for metastatic castration-resistant prostate cancer: a cohort study within the French nationwide claims database.

OBJECTIVES: To conduct the direct comparison of abiraterone acetate and docetaxel for first-line treatment of metastatic castration-resistant prostate cancer (mCRPC) in real-life settings. METHODS: Data were extracted from the French nationwide claims database (SNDS) on all men aged ≥40 years starting first-line treatment with abiraterone acetate or docetaxel for mCRPC in 2014. A high-dimensional propensity score including 100 baseline characteristics was used to match patients of both groups and form two comparative cohorts. Three-year overall survival and treatment discontinuation-free survival were determined using Kaplan-Meier analysis. RESULTS: In 2014, 2,444 patients started abiraterone for treatment of mCRPC and 1,214 started docetaxel. After trimming and matching, 716 patients were available in each group. Median overall survival tended to be longer in the abiraterone acetate cohort (23.8 months, 95% confidence interval = [21.5; 26.0]) than in the docetaxel cohort (20.3 [18.4; 21.6] months). Survival at 36 months was 34.6% for abiraterone acetate and 27.9% for docetaxel (p = 0.0027). Treatment discontinuation-free median was longer in the abiraterone acetate cohort compared to the docetaxel cohort (10.8 [10.1; 11.7] versus 7.4 [7.0; 8.0] months). CONCLUSION: The findings underline the interest of oral abiraterone acetate over intravenous docetaxel as the first-line treatment option in mCRPC.

G-Quadruplex Recognition by DARPIns through Epitope/Paratope Analogy.

We investigated the mechanisms leading to the specific recognition of Guanine Guadruplex (G4) by DARPins peptides, which can lead to the design of G4 s specific sensors. To this end we carried out all-atom molecular dynamic simulations to unravel the interactions between specific nucleic acids, including human-telomeric (h-telo), Bcl-2, and c-Myc, with different peptides, forming a DARPin/G4 complex. By comparing the sequences of DARPin with that of a peptide known for its high affinity for c-Myc, we show that the recognition cannot be ascribed to sequence similarity but, instead, depends on the complementarity between the three-dimensional arrangement of the molecular fragments involved: the α-helix/loops domain of DARPin and the G4 backbone. Our results reveal that DARPins tertiary structure presents a charged hollow region in which G4 can be hosted, thus the more complementary the structural shapes, the more stable the interaction.

How Fragile We Are: Influence of Stimulator of Interferon Genes (STING) Variants on Pathogen Recognition and Immune Response Efficiency.

The stimulator of interferon genes (STING) protein is a cornerstone of the human immune response. Its activation by cGAMP in the presence of cytosolic DNA stimulates the production of type I interferons and inflammatory cytokines. In the human population, several STING variants exist and exhibit dramatic differences in their activity, impacting the efficiency of the host defense against infections. Understanding the molecular mechanisms of these variants opens perspectives for personalized medicine treatments against diseases such as viral infections, cancers, or autoinflammatory diseases. Through microsecond-scale molecular modeling simulations, contact analyses, and machine learning techniques, we reveal the dynamic behavior of four STING variants (wild type, G230A, R293Q, and G230A/R293Q) and rationalize the variability of efficiency observed experimentally. Our results show that the decrease in STING activity is linked to a stiffening of key structural elements of the binding cavity together with changes in the interaction patterns within the protein.

Specific Recognition of the 5'-Untranslated Region of West Nile Virus Genome by Human Innate Immune System.

In the last few years, the sudden outbreak of COVID-19 caused by SARS-CoV-2 proved the crucial importance of understanding how emerging viruses work and proliferate, in order to avoid the repetition of such a dramatic sanitary situation with unprecedented social and economic costs. West Nile Virus is a mosquito-borne pathogen that can spread to humans and induce severe neurological problems. This RNA virus caused recent remarkable outbreaks, notably in Europe, highlighting the need to investigate the molecular mechanisms of its infection process in order to design and propose efficient antivirals. Here, we resort to all-atom Molecular Dynamics simulations to characterize the structure of the 5'-untranslated region of the West Nile Virus genome and its specific recognition by the human innate immune system via oligoadenylate synthetase. Our simulations allowed us to map the interaction network between the viral RNA and the host protein, which drives its specific recognition and triggers the host immune response. These results may provide fundamental knowledge that can assist further antivirals' design, including therapeutic RNA strategies.

Autophagy and evasion of the immune system by SARS-CoV-2. Structural features of the non-structural protein 6 from wild type and Omicron viral strains interacting with a model lipid bilayer.

The viral cycle of SARS-CoV-2 is based on a complex interplay with the cellular machinery, which is mediated by specific proteins eluding or hijacking the cellular defense mechanisms. Among the complex pathways induced by the viral infection, autophagy is particularly crucial and is strongly influenced by the action of the non-structural protein 6 (Nsp6) interacting with the endoplasmic reticulum membrane. Importantly, differently from other non-structural proteins, Nsp6 is mutated in the recently emerged Omicron variant, suggesting a possible different role of autophagy. In this contribution we explore, for the first time, the structural properties of Nsp6 thanks to long-timescale molecular dynamics simulations and machine learning analysis, identifying the interaction patterns with the lipid membrane. We also show how the mutation brought by the Omicron variant may indeed modify some of the specific interactions, and more particularly help anchor the viral protein to the lipid bilayer interface.

Never Cared for What They Do: High Structural Stability of Guanine-Quadruplexes in the Presence of Strand-Break Damage.

DNA integrity is an important factor that assures genome stability and, more generally, the viability of cells and organisms. In the presence of DNA damage, the normal cell cycle is perturbed when cells activate their repair processes. Although efficient, the repair system is not always able to ensure complete restoration of gene integrity. In these cases, mutations not only may occur, but the accumulation of lesions can either lead to carcinogenesis or reach a threshold that induces apoptosis and programmed cell death. Among the different types of DNA lesions, strand breaks produced by ionizing radiation are the most toxic due to the inherent difficultly of repair, which may lead to genomic instability. In this article we show, by using classical molecular simulation techniques, that compared to canonical double-helical B-DNA, guanine-quadruplex (G4) arrangements show remarkable structural stability, even in the presence of two strand breaks. Since G4-DNA is recognized for its regulatory roles in cell senescence and gene expression, including oncogenes, this stability may be related to an evolutionary cellular response aimed at minimizing the effects of ionizing radiation.